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We studied the photosynthetic activity of Scots pine (Pinus sylvestris L.) and Norway 
spruce (Picea abies [L] Karst) in relation to air temperature changes from March 2013 
to February 2014. We measured the chlorophyll fluorescence of approximately 50 trees 
of each species growing in southern Finland. Fluorescence was measured 1-3 times 
per week. We began by measuring shoots present in late winter (i.e., March 2013) 
before including new shoots once they started to elongate in spring. By July, when the 
spring shoots had achieved similar fluorescence levels to the older ones, we proceeded 
to measure the new shoots only. We analyzed the data by fitting a sigmoidal model 
containing four parameters to link sliding averages of temperature and fluorescence. A 
parameter defining the temperature range over which predicted fluorescence increased 
most rapidly was the most informative with in describing temperature dependence of 
fluorescence. The model generated similar fluorescence patterns for both species, but 
differences were observed for critical temperature and needle age. Down regulation of 
the light reaction was stronger in spring than in autumn. Pine showed more conservative 
control of the photosynthetic light reactions, which were activated later in spring and more 
readily attenuated in autumn. Under the assumption of a close correlation of fluorescence 
and photosynthesis, spruce should therefore benefit more than pine from the increased 
photosynthetic potential during warmer springs, but be more likely to suffer frost damage 
with a sudden cooling following a warm period. The winter of 2013-2014 was unusually 
mild and similar to future conditions predicted by global climate models. During the mild 
winter, the activity of photosynthetic light reactions of both conifers, especially spruce, 
remained high. Because light levels during winter are too low for photosynthesis, this 
activity may translate to a net carbon loss due to respiration. 
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INTRODUCTION 

Trees growing in the boreal zone must respond to a wide range 
of climatic conditions during a typical year. For evergreens, this 
means that photosynthetic organs must be active in the sum- 
mer but become dormant and frost-resistant for the winter 
(Leinonen and Hanninen, 2002). Spring is a demanding time for 
boreal conifers when light is abundant but cold air temperatures 
and possibly lack of available water can hinder photosynthesis 
(Suni et al., 2003). Trees are locally adapted to prevailing con- 
ditions, but global warming will change the synchronization of 
temperature and light regime at a given location, i.e., although 
temperatures are expected to be higher, incident radiation will 
remain the same. Increasing temperatures will affect respira- 
tion during autumn and winter (Vesala et al., 2010) and pos- 
sibly enhance photosynthetic activity during spring (Hall et al., 
2009). Therefore, climate change will impact the productivity of 
boreal forests (Wang et al., 2011), but the nature of that impact 



will depend on how trees adapt to the new temperature-light 
regime (TLR). 

The photosynthetic pathway in coniferous trees is typical and 
consists of two parts; a light reaction captures photons and stores 
their energy in short-lived compounds before a dark (i.e., light 
independent) reaction converts them into more stable products. 
The efficiency of the light reaction is largely determined by inci- 
dent radiation while the dark reaction is a complex biochemical 
process more dependent on prevailing temperatures (Porcar- 
Castell et al, 2008). Consequently, a cold but bright TLR could 
lead to unchecked light reactions producing free radicals that 
damage plant tissues (Huner et al., 1998; Porcar-Castell et al., 
2008). In response to this problem, boreal evergreens can actively 
downgrade light reactions taking place in their photosynthetic 
tissues (Kolari et al, 2007). 

As expected for locally-adapted species, the potential photo- 
synthesis of evergreens follows local changes in the TLR (Ottander 
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and Oquist, 1991; Beck et al., 2004). Due to the asynchrony 
of temperature and light, the phenology of a given species will 
affect its photosynthesis, growth, reproduction, and risk of frost 
damage in altered climates of the future. Development of the 
photosynthetic capacity of boreal conifer species has been exten- 
sively studied (Bergh et al., 1998; Lundmark et al., 1998; Leinonen 
and Hanninen, 2002; Makela et al, 2004; Kolari et al, 2007), 
but comparative studies with different species in the same envi- 
ronmental conditions are limited (Lundmark et al., 1988; Ogren 
et al., 1997). Clearly, if conifers show species-specific responses 
to a given change in TLR, consequent changes to the produc- 
tivity of boreal forests in the future will be determined by the 
co-distribution of species and the new climate. Comparative anal- 
yses of species living in similar conditions and modeling the data 
appropriately would improve the precision with which impacts 
of climate change on productivity and ecosystem health can be 
predicted. However, observation periods must include a wide 
range of climatic conditions in order to detect and investigate any 
differences among species. 

The goal of this study was to understand the light reaction 
dynamics of Norway spruce and Scots pine in response to tem- 
perature variation during the year. We also compared chlorophyll 
fluorescence between spring and autumn. We hypothesized that 
in comparison to Scots pine, photosynthetic light reactions in 
Norway spruce resume earlier in spring and remain at higher 
levels for most of the winter. We fitted nonlinear mixed-effect 
models to light reaction activity based on average and minimum 
daily temperatures, and compared the predictions for spruce and 
pine in relation to changes in TLR. 

MATERIALS AND METHODS 

We measured the development of fluorescence in 51 Norway 
spruce (Picea abies [L.] Karst) and 48 Scots pine (Pinus sylvestris 
L.) located in the Haapastensyrja common garden (Finnish Forest 
Research Institute) in southern Finland (60°37', 24°25', 1250 
d.d.). Most sample trees originate from Finland (60-70°N) but 
some of the spruce trees came from populations in central Europe 
(Table 1). Norway spruce are growing in a stand planted in 
1980 and were approximately 15 m tall when measurements were 
taken. For practical reasons, data were collected from shoots that 
were on the outer edges of the spruce stand where trees had 
healthy shoots at least down to breast-height level. We assumed 
this to be an indication that these shoots received sufficient light 
to maintain good condition. 

The Scots pine were grafted from shoots collected from mature 
trees in spring 2006 to rootstocks sown a year before. The grafts 
were planted in an open field adjacent to the spruce stand in 2007. 
In 2013, the trees were about 1.5-2.0 m tall; data were collected 
from the upper branches. All Scots pine originated from northern 
Finland (Table 1). 

Fluorescence measurements began 06 March 2013 and con- 
tinued until the end of February 2014. Measurements were 
generally made weekly, increasing to 2-3 times a week when 
abrupt weather changes and consecutive changes in the fluores- 
cent activity of the trees were expected. Measurements began with 
needles formed during 2012 (hereafter referred to as 2012 nee- 
dles), and were later expanded to include needles that formed 



Table 1 | Origin and number of trees studied in this experiment. 





Origin 


Number of trees 


Latitude, °IM 




Switzerland 


1 


46.8 




Czechoslovakia 


4 


48.9 


CO 


Poland 


2 


51 .4 


CD 


Loppi, Finland 


4 


60.7 


Norway spruce 
(Picea abies [L] 


Langelmaki, Finland 
Muurame, Finland 
Nurmes, Finland 
Puolanka, Finland 
Simo, Finland 


5 
6 
5 
5 
11 


61 .6 
61.9 
63.5 
64.9 
65.7 




Ylitornio, Finland 


4 


66.3 




Sodankyla, Finland 


2 


67.4 




Muonio, Finland 


2 


68.0 




Total 


51 













_i 








■23 


Ivalo, Finland 


13 


68.6 




Kevo, Finland 


12 


69.8 


CD -5 






c ^ 
•q. °> 


Rovaniemi, Finland 


11 


66.4 


ots 
inus 


Sodankyla, Finland 


12 


67.3 




Total 


48 





in the early summer of 2013 (hereafter 2013 needles), once they 
had bursted. Both 2012 and 2013 needles were measured in par- 
allel until their fluorescence was the same in early luly, after 
which only 2013 needles were measured. Measuring needles of 
different ages enabled the development of fluorescence between 
spring (i.e., 2012 needles) and autumn (i.e., 2013 needles) to be 
observed. Fluorescence was measured in the field with a hand- 
held Pocket PEA fiuorimeter (Hansatech, Ltd.), where needles 
were first enclosed in light-proof clips for at least 30 mins to allow 
all photosynthetic reactions time to complete. The fiuorimeter 
was then attached to the clip and a strong pulse of light was passed 
through the sample. Initially, most of the light energy is absorbed 
and processed by the dark reaction and the fluorescent irra- 
diance is at its lowest value (F 0 ). Shortly thereafter the dark 
reaction becomes saturated and surplus energy from the light 
reaction causes the fluorescent irradiance to reach a peak value 
(-Fm). The fiuorimeter records the initial and saturated fluorescent 
irradiance and calculates a relative index as: 



Fo 



Fv 



(1) 



The expression F v /F m reflects the photosynthetic potential: when 
light reactions are suspended, there is minimal change in the flu- 
orescent irradiance and, therefore, values are small, while high 
values indicate that the photosynthetic apparatus is active. 

Leaf clips for the Pocket PEA fiuorimeter are intended to mea- 
sure broadleaf species with flat leaves. We modified the clips 
to ensure a light-proof closure of the clip around the cylin- 
drical sample needles (Figure 1). For each tree, one shoot was 
selected and tagged for measurement. Typically, two needles were 
inserted into the aperture of the clip. In repeated measurements, 
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placement of the clip on the shoot was selected randomly and 
independently of earlier samples. 

Throughout the measuring period, outdoor temperatures were 
recorded with a NTC-type thermistor located in a weather-logger 
box in an open field about 200 m from the stands. Temperature 
was measured once per minute. Daily average temperature was 
calculated as the arithmetic average of all daily measurements, 




FIGURE 1 | Clips used to measure chlorophyll fluorescence modified 
with an extra foam pad added to the lower jaw and a small rubber 
band to draw the clip closed. These modifications helped provide a 
light-proof closure of the clip around the needles to be measured. 



and daily minimum temperature was produced by selecting the 
minimum value for each day (Figure 2). To account for the time 
lag in the reaction of trees to changing conditions, we calcu- 
lated a moving average of the temperature using a range of time 
constants from 1 to 12 days based on values reported in the 
literature (see discussion). For comparison (Figure 2), we calcu- 
lated a long-term daily average for 1961-2012 from temperature 
measurements provided by the Finnish Meteorological Institute 
(Venalainen et al., 2005). 

The dependence of fluorescence measurements (photosyn- 
thetic activity) on temperature was described with a non-linear 
mixed effects model based on a four-parameter logistic response 
function: 

~ , #max #min 

F=a min -\ jr- =-, (2) 

1 + exp[ 

where F is the predicted fluorescence, T is the sliding average 
of either the daily average or daily minimum temperature,a m ,„ 
and a max are the minimum and maximum asymptotes, respec- 
tively, b is a scale parameter determining the rate of ascent (i.e., 
a lower value of b translates to a steeper rate of ascent), and c 
is the temperature at which F is halfway between a m \ n and a max - 
This temperature value also corresponds to the largest change of 
fluorescence activity. Variable c is referred to as the crossover tem- 
perature, interpreted as the temperature value characterizing the 
critical temperature range for photosynthetic light reactions. The 
asymptotic values a mm and a max were fitted as fixed parameters 



30 




1-Jan 1-Feb 1-Mar 1-Apr 1-May 1-Jun 1-Jul 1-Aug 1-Sep 1-Oct 1-Nov 1-Dec 1-Jan 1-Feb 

Date 



FIGURE 2 1 Average (red line) and minimum (blue line) daily 
temperatures during the study period. Stronger colors indicate 
the period of fluorescence measurements, with 2 months 
preceding the period presented in shaded colors. The thinner red 



and blue lines show the long-term mean and minimum daily 
temperatures (red and blue, respectively) from the Finnish 
Meteorological Institute weather data from the years 1961 to 2012 
(Venalainen et al., 2005). 
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separately for each species. The other two parameters, b and c, 
were fitted separately for each species as well as for the set of 
observations representing spring and autumn (i.e., 2012 and 2013 
needles ) . Furthermore, b and c were allowed to vary between trees 
so that the value of b for the ith observation is as follows: 



h = &s(j), yd) + Ph 



(3) 



where s(i) and y(i) indicate the species and needle cohort, respec- 
tively, and p\ is a zero-mean tree-specific random effect; c, was 
modeled in a similar fashion. Variances for the random effect were 
estimated separately for each species, but a common estimate for 
residual variance was appropriate. Based on the latter assumption, 
both species could be included in the same model, allowing for 
the statistical testing of differences between them. Model fitting 
was carried out in the R-environment (R Core Team, 2013) with 
the nlme package (Pinheiro et al, 2013), which uses approximate 
maximum likelihood (Lindstrom and Bates, 1990). For Norway 
spruce, association between tree origin and cross-over temper- 
ature was assessed with Spearman's rank correlation (Best and 
Roberts, 1975) of latitude of origin and estimated tree-specific 
random effect in c;. 

RESULTS 

Fluorescence values (Figure 3) showed that in spring 2013 the 
activity of the photosynthetic light reactions in Norway spruce 
started to develop earlier, (i.e., at lower air temperatures) than 
in Scots pine. Fluorescence in Scots pine was lower than that of 



spruce throughout the 6-week period of development of photo- 
synthesis in the spring. Both pine reached the same fluorescence 
level as spruce by the end of May. 

During the autumn, Scots pine began attenuating light reac- 
tions earlier and to lower levels compared to Norway spruce, 
but as long as autumn temperatures remained mild, fluores- 
cence activity did not reach the lowest levels observed in the 
previous spring. The winter of 2013-2014 was unusually mild in 
southern Finland, with only mildly freezing temperatures except 
for a cold period of 3 days in early December and another 
of 2 weeks in late January (Figure 1). During the first cold 
spell, pine fluorescence decreased to a much lower level than 
spruce. In both species, the photosynthetic apparatus was reac- 
tivated after the cold spell, but pine resumed at a lower level 
than spruce. The cold period in January 2014 caused pine flu- 
orescence to decrease to a similar level to the lowest observed 
in early spring 2013, while spruce persisted at a higher level. 
Light reactions of both species reactivated when air temperatures 
warmed, but pine resumed activity at a lower level compared to 
spruce. 

A range of time constant values for the moving average of 
temperature was tested, and values were compared by fitting the 
mixed model to each temperature record. The log-likelihood val- 
ues of the model for each value of the time constant are given 
in Table 2. A time constant value of 7 days gave the largest log- 
likelihood, but differences among adjacent values were small. 
Also, daily mean or daily minimum temperature only had a slight 
effect on fit of the mixed model (Table 2). 



1 

0.9 
0.8 



_ 0.7 

E 



!±- 0.6 



C 

8 0.5 



o 

3 0.4 



0.3 
0.2 
0.1 
0 




re 
2 




2 

i 

r-i 

ro 



O 



-o- Spruce 2012 needles 
— o - Spruce 2013 needles 
-o Pine 2012 needles 
---■o- Pine 2013 needles 



| 

00 



4* 



3 



Date 



FIGURE 3 | Development of chlorophyll fluorescence March 
2013-February 2014. Darker blue and red lines with circles indicate 
measurements for needles that burst in spring 2012 for pine and spruce, 



respectively. Brighter lines with diamonds are the corresponding values for 
needles that burst in spring 2013. Error bars show the standard deviation of 
measurements for each day. 
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The crossover temperature parameter value, which indicates 
the temperature where the change of fluorescent activity is at its 
largest, was lower for spruce compared to pine and also lower 
for spring needles (or "2012 needles") compared to autumn nee- 
dles ("2013 needles") (Table 3, Figure 4). This means that for any 
given temperature below the saturating temperature, the photo- 
synthetic activity of spruce is higher than that of pine. Needles 
were more active in autumn than during spring for both species. 
The slope of fluorescence development was steeper for pine in the 
spring. This suggests a more rapid development of photosynthetic 
capacity as temperatures increase, which compensates for some 
of the effect of the higher crossover temperature. According to 
the fitted models, differences in all other parameter values except 
for a m i n between species and needle origins are clearly statistically 
significant (p-values < 0.0001). 

Models for individual trees show greater variation for spruce 
than pine (Table 3, Figure 4). Spruce trees were more variable 
in terms of their origin, being drawn from throughout Finland 
(i.e., 60-66°N) and locations in central Europe, while pines 
all came from northern Finland (Table 1). However, association 



Table 2 | The log-likelihood values of the mixed model with different 
time constant values for the moving mean and minimum daily 
temperatures. 



Time constant Log-likelihood 



Days 


Mean temperature 


Minimum temperature 


1 


12.771 


9428 


2 


14.029 


11.476 


3 


14.359 


12.546 


4 


14.545 


13.483 


5 


14.648 


14.094 


6 


14.692 


14.339 


7 


14.695 


14.432 


8 


14.669 


14.457 


9 


14.621 


14.439 


10 


14.559 


14.393 


11 


14.486 


14.326 


12 


14.407 


14.243 



between the latitude of origin and cross-over temperature was not 
statistically significant in our data (Figure 5). 

DISCUSSION 

Our measurements of chlorophyll fluorescence indicate that the 
activity of photosynthetic light reactions of boreal conifers fol- 
lows ambient air temperature. In spring 2013, F v /F m reached 






Temperature (°C) 

FIGURE 4 | Predicted chlorophyll fluorescence (F„/F m ) as a function of 
temperature for the spruce and pine needles (upper and lower rows, 
respectively) and for the old and new needles (left and right columns, 
respectively). The thick lines show the average model for all trees in the 
group, with thin lines indicating individual trees. Figures are based on the 
best-fitting model using a 7-day sliding average of daily mean temperature. 



O 

o 

CD 

D 

CO 
s— 
CD 
CL 



CD 
> 
O 
in 

C/) 

2 
o 



CM- 



— i 1 1 r~ 

50 55 60 65 
Latitude (°N) 



FIGURE 5 | Association between the latitude of origin and estimated 
tree-specific cross-over temperature c, for spruces in the spring (2012 
needles). Spearman's rank correlation p = —0.23 (p= 0.10). 



Table 3 | Parameter estimates and standard errors (in parentheses) 
for the fixed effects in the non-linear mixed model using a 7-day 
sliding average of daily mean temperature as a predictor. 

Pine [Pinus Spruce (Picea abies 

sylvestris L) [L] Karst) 

a max 0.847(0.002) 0.833(0.001) 

a min 0.081 (0.003) 0.078 (0.006) 

cm.1 c 3.841(0.081) -0.925(0.318) 

§c§" b 2.181(0.062) 3.834(0.095) 



c 

| c -1.647(0.119) -10.765(0.341) 

5 1 b 4.886(0.106) 4.710(0.122) 

CM IS 
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minimum levels in March after a cold period of several weeks. 
Air temperatures late in 2013 and during January-February 2014 
were unusually warm, up to 6°C above the long-term aver- 
age (Figure 1). During this period, F v IF m values remained high, 
with spruce almost reaching levels observed during the sum- 
mer (Figure 3). Studies of net ecosystem exchange using eddy- 
covariance methods (Van Dijk and Dolman, 2004; Vesala et al., 
2010) have shown that trees maintain their respiration during 
such warm but dark periods, which suggests that photosynthetic 
dark reactions are also active during warm spells in winter. The 
low F v IF m values observed after the cold period in January 2014 
indicate a down-regulation of light reactions. Orlander (1993) 
suggested that this down-regulation is due to damage to photo- 
system II caused by light reception in subfreezing temperatures. 
However, the Norway spruce and Scots pine in our common 
garden experiment quickly recovered from this decline, and the 
F v IF m ratio followed the temperature increase. This indicates that 
photosystem II of the conifers can sustain freeze-thaw cycles of at 
least -20° C, with rapid recovery during milder periods. Norway 
spruce tended to maintain a higher F v /F m ratio, and reactivated 
as rapidly as Scots pine, which showed stronger photoinhibition 
during cold spells. Differences in the activity of the photosyn- 
thetic apparatus can be linked to the way these two conifers 
down-regulate respiration. Ogren et al. (1997) showed that cold 
tolerance of conifers is correlated with the soluble sugar content 
of their needles and that Norway spruce is more frost tolerant 



than Scots pine or Lodgepole pine — a proposal that is in conflict 
with results presented here. The authors also observed that spruce 
down-regulated respiration faster and more extensively than pine, 
which improved the conservation of sugar storage in spruce. 

During the mild winter period, our spruce trees exhibited con- 
siderably higher values of fluorescence compared to pine. These 
results reflect those of Beuker et al. (1998), who tested frost tol- 
erance by exposing seedlings to a range of freezing temperatures 
and measuring the subsequent leakage of electrolytes. They found 
that pine becomes frost tolerant much earlier during the autumn, 
loses it later in the spring, and tolerates lower temperatures dur- 
ing the winter (Figure 6, a reproduction of Figures 3, 6 in Beuker 
et al., 1998). Several experimental studies have shown that frost 
resistance and photosynthetic activity are linked (Leinonen, 1996; 
Repo et al., 2000, 2006; Makela et al, 2004). Therefore, the con- 
clusions of Beuker et al. (1998) are consistent with our results. 
We did not measure actual needle temperatures, as Martin et al. 
(1999) showed that needle temperatures of subalpine Pacific sil- 
ver fir are typically < 1°C higher than ambient, and we assumed a 
similar relationship was applicable in southern Finland. 

In this study, the crossover parameter values (Table 3, 
Figure 3) were lower in the spring (i.e., 2012 needles) compared 
to the autumn (i.e., 2013 needles) for both species. This indicates 
a lower activity of light reactions at any temperature in the spring 
compared to the autumn. This agrees with Repo et al. (2006), 
who noted that the temperature-linked photosynthetic activity is 
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FIGURE 6 | Development of frost resistance of Norway spruce (red) and 
Scots pine (blue) buds through three winters. The solid line shows the 
daily mean temperature for the study period. The figure is a combination of 



Figures 3, 6 in Beuker et al. (1998) (reproduced with permission), who tested 
frost resistance by exposing seedlings to freezing temperatures and 
measuring the subsequent leakage of electrolytes. 
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different in autumn and spring and that photosynthetic poten- 
tial is higher in autumn compared to spring, when estimated with 
light-saturated assimilation rate and the apparent quantum yield 
achieved with gas exchange measurements. 

Results presented here indicate that Norway spruce regains 
the potential for photosynthetic activity more rapidly than Scots 
pine when temperatures suddenly warm during the cold months. 
The results disagree with those of Lundmark et al. (1988), who 
suggested these species to be similar in this respect. Norway 
spruce is considered to be more prone to spring frost damage 
due to a preference for frost-prone sites and sensitivity of new 
shoots (Lundmark and Hallgren, 1987). Our results indicate that 
this sensitivity is also a consequence of earlier reactivation of 
photosynthesis. 

Because the Scots pine in our study had a more northerly ori- 
gin than the Norway spruce, we expected the pines to be better 
adapted to cold and that they would start photosynthesizing at 
lower temperatures (Hanninen et al., 2007). However, our results 
indicate that photosynthetic reactivation was actually more con- 
servative in pines. We conclude that the difference between these 
two conifers is due to species-specific responses to a common 
environmental change rather than a consequence of sample tree 
provenance. Moreover, spruce trees were drawn from a wider area 
than the pines, so one would assume origin explains some of the 
variance in the model parameters for this species. However, we 
found no statistically significant relationship between tree origin 
and parameter values, and an explanation of the greater variation 
observed for spruce remains obscure. 

We tested a range of time constant values from 1 to 12 days 
when calculating the moving average of temperature. In the end, 
the time constant had only a slight effect on model performance. 
The time constant, or response time, is believed to reflect the rate 
at which photosynthetic activity reacts to changes in the envi- 
ronment. Earlier studies have produced a wide range of response 
times, e.g., 2 days (Ottander and Oquist, 1991), 8 days (Kolari 
et al., 2007), or up to 12 days (Makela et al., 2004). Our model fit- 
ted the data slightly better with a response time of 7 days, but with 
only slight differences in model performance between one value 
and another, it is difficult to select any specific value. We suspect 
that the flat distribution of model performance over a range of 
response times accounts for the wide range of values reported in 
the literature. 

Differences in the phenology of photosynthesis between the 
two conifer species have important implications for modeling for- 
est productivity and forestry planning in the future. According 
to our results, spruce reactivates the light reactions of its photo- 
synthesis more easily. We speculate that this implies that spruce 
also more readily reactivates its photosynthesis, and is therefore 
better adapted to utilize the warmer spring temperatures. On the 
other hand, earlier onset of photosynthetic activity also runs a 
greater risk of frost damage, especially if variation in springtime 
daily temperatures increases. However, if daily temperature vari- 
ation remains similar to that seen today, Norway spruce will gain 
an advantage in a warmer future climate. Our model of photosyn- 
thetic activity can be used in process-based growth/yield models 
to characterize the different patterns of photosynthesis and res- 
piration for spruce and pine. It can also be used as a tool for 



exploring how the risk of frost will develop under future climate 
scenarios and weather patterns. 
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